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Abstract 

We study pion production from proton synchrotron radiation in the presence of strong magnetic 
fields. We derive the exact proton propagator from the Dirac equation in a strong magnetic field 
by explicitly including the anomalous magnetic moment. In this exact quantum-field approach the 
magnitude of pion synchrotron emission turns out to be much smaller than that obtained in the 
semi-classical approach. However, we also find that the anomalous magnetic moment of the proton 
greatly enhances the production rate about by two order magnitude. 

PACS numbers: 95.85.Ry,24.10.Jv,97.60.Jd, 
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I. INTRODUCTION 


Magnetic fields in nentron stars play an important role in the interpretation of many 
observed phenomena. Indeed, strongly magnetized neutron stars (dubbed magnetars QS) 
hold the key to understanding the asymmetry in supernova (SN) remnants and the still 
unresolved mechanism for non-spherical SN explosions. Such strong magnetic helds are also 
closely related to the unknown origin of the kick velocity that proto-neutron stars (PNSs) 
receive at birth. 

It is widely accepted that soft gamma repeaters (SGRs) and anomalous X-ray pulsars 
(AXPs) correspond to magnetars [^, and that the associate strong magnetic helds have a 
signihcant role in production high energy photons. Furthermore, short duration gamma-ray 
bursts (GRB^ may arise from highly magnetized neutron stars 5| or mergers of binary neu¬ 
tron stars 6|-[8|, and the most popular theoretical models for the long-duration GRBs 9l-ll2| 
invoke magnetized accretion disks around neutron stars or rotating black holes (collapsars) 
for their central engines. Such magnetars (or black holes with strong magnetic helds) have 
also been proposed 3, 141 as an acceleration site for ultra high-energy (UHE) cosmic rays 
(UHECRs) and a possible association between magnetar hares 1^ and UHECRs has also 
been observed. 

When a particle is accelerated in an external held or in collision with another particle, it 
can emit quanta corresponding to the held with which the particle interacts. Synchrotron 
radiation can be produced by high-energy protons accelerated in an environment containing 
a strong magnetic held. This process has been proposed as a source for high-energy photons 


in the GeV — TeV range 


16 


2l[ |. possibly in association with GRBs. Synchrotron emission. 


however, can occur through any quanta that may couple to an accelerated particle. Since 
protons strongly couple to meson helds, a high-energy proton can also radiate pions and 
other mesons, as well as photons. 

In fact, the meson-nucleon couplings are about 100 times larger than the photon-nucleon 
coupling, and the meson production process is expected to exceed photon synchrotron emis¬ 


sion in the high energy regime. For example. Refs. 


22 


26| addressed the possibility of 7r° 


emission from a proton in a strong magnetic held. However, these calculations were per¬ 
formed approximately; i.e. in a semi-classical way and an approximate quantum-mechanical 
treatment of the proton transitions among the Landau levels of the strong magnetic held. 
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In the semi-classical approach for the production of synchrotron radiation, the magnetic 
field strength is characterized by the curvature parameter, y = Ef/{m^Rc), given in terms of 
the incident particle mass m, its energy Ei, and the curvature radius Re- For y << 1 (i?c >> 
1), the radiation can be treated in a classical way. For example, in a gravitation field or 
a relatively small electro-magnetic held, charged particles may have such large curvature 
radii. However, for y >> 1 (i?c << 1) he., very high energetic particles propagating in 
the very small curvature radius of a strong magnetic held, quantum ehects must be taken 
into account. Since we consider a proton propagating in a strong magnetic held, which may 
imply a small Lamor radius Ri = E/{eB), the radiation should be also treated in a quantum 
mechanical way along with the Landau quantization by the magnetic held. 

For proton propagation in a strong magnetic held, the value of y can be given as 


where Bcr = M^/e and B± is the magnetic held perpendicular to the proton momentum 
direction. For illustration, we assume a 1 ~ 5 GeV proton propagating in a strong magnetic 
held with B ~ 10^®G, for which y ~ 0.01. Through analysis of the pion emission based on a 
quantum held theory treatment, we show that the pion emission can be comparable to the 
photon emission. In particular, we demonstrate that the anomalous magnetic moment of 
the charged particle plays a vital roles in the pion synchrotron emission. 



Ef _eE,B_E,B^ 
M^Rl M3 Mp^B,r 


In this work, we exploit the Green’s function method for the propagation of protons 
in a strong magnetic held. This approach has not previously been applied to derive pion 
synchrotron emission. For the pion-nucleon coupling, the p-wave interaction is dominant. 


and the pion emission amp 
emitted pion momentum Q 


itude is mainly proportional to < 'ijjj^cr ■q'lpN >, with q being the 
Because the pion emission along the direction of the magnetic 
held is not allowed by the conservation of energy-momentum, < 'ip^(j±'ipN > has a dominant 
contribution. Hence, the spin-hip contributions necessary for pseudo-scalar emission are 
signiheant. 

Section 2 is devoted to an introduction of our theoretical formalism based upon the 
Green’s function method. Numerical results are presented in section 3 along with detailed 
discussions. A summary and conclusions are presented in section 4. 


^ Note: this fact is independent of the choice of whether one has PS- or PV-coupling. 
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II. FORMALISM 


A. Proton Green’s Pnnction 

We assume a uniform magnetic field along the 2 ;-direction as B = (0, 0, B) and take the 
electro-magnetic vector potential to be A = (0, 0, xB, 0) at the position r = {x, y, z) . 
The relativistic proton wave function 'ijj is obtained from the following Dirac equation; 


jf, ■ {id^ - eA^) - tun 


eKr, 


2m AT 




'^{x) = 0, 


( 2 ) 


where m^r is the proton mass, Kp is the proton anomalous magnetic moment (AMM), and e 
is the elementary charge. 

Here, we scale all variables with \/eH as and Mjv = and we 

write the wave function, e.g. 


^(X) = 


J{PyX2+P.X3-iEXo) 


( \rUi{X-Py)\ 

As UX-Py) 

X^fn+liX-Py) 

\X, UX-Py) J 

Eq. ([2]) then leads the following characteristic equation; 

E70 ^J2{n + l)'y‘^ - PU - {hp/Mn)E^ - 


( 3 ) 


A = 0 


( 4 ) 


with A = (Ai, As, A3, A4) and = diag{l, —1,1, —1). 

Here, the Landau level n starts from n = —1 when s = 1 and from n = 0 when s = —1. 
We then redehne the Landau number ul s.s hl = n + [1 + s)/2. By solving Eq. (jlj), we then 
obtain the energy eigenvalues as 


EijiL, Pz, s) = ±y (^^2nL + l- s + Mjf - sKp/MU. 

The proton Green’s Function G in a magnetic held is written as 

Pm\^l,s,Pz) , p^^\nL.x.P, 


( 5 ) 


ni^=0 s=±l 


-7 + 


Po - E{nL, s, P^) +iS Po + E{nL, s, P^) id 


F{X') (6) 


with 


F = - P,)pp + (A- - 


(7) 


4 


















( 8 ) 


^ [E'jo + \/2nL + 1 - S7^ - + Mn + {hp/Mn)T.^] 


X 1 + 


stip/Mj 


N 


\j2nL + 1 — s + 


+ 75^ L 


Pm = ^ [^7o - \/2ni, + 1 - S7^ + P 27 ^ - - {kp/Mn)E^] 

SKp/M^ 


X 1 + 


- 75^ , 


^j2,nL + 1 — s + 

where Ep is the Fermi energy, and the spin vector a is dehned as 

sPr „ „ sE 


a(p, s) = 


= , 0 , 0 , 


\j2np + 1 — s + Mfq j 


(9) 


( 10 ) 


B. Pion Production 


In this subsection we consider the pion production rate in the presence of a strong mag¬ 
netic held. We start the calculation from the following pseudo-vector coupling interaction 
Lagrangian density as 

C = —'ll!-f5-f^Ta'lljd^4>a, (H) 

where is the pseudo-vector pion-nucleon coupling constant, is the pion mass, and 
(j) is the pion held. We then calculate the pion decay rate from an initial proton with 
np = rii, s = Si and Pz = Piz to a hnal proton with np = n/, s = Sf and Pz = Pfz- The 
pion momentum scaled by is written as Q = (P^r, 0, Qr, Qz), where without loss of 

generality the transverse pion momentum is assumed to be directed along the y-a.xis. 

Using the proton propagator, the proton self-energy from one-pion exchange can be writ¬ 
ten as ^ 

S.(Pl, P 2 ) = t WTaGiRl, P2)ra7.755^-9"A^(Pi - P 2 ), (12) 

and the decay width is calculated from the imaginary part of the expectation value of the 
self-energy as 

=-Im J d^Rid'^R 2 i>i{Ri)T.^{Ri,R 2 )i>iiR 2 ), (13) 


where is the pion propagator. 

By performing a Fourier transformation, we can obtain the diherential decay width of 
the proton as 


d^TpjV^ ^ 1 ^ 6{Ef + E^-E,) 

dQ^ Stt^E^ [m^J ^ AEiEf "R 

^ ' rif.Sf 


(14) 
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with 

Wif = AEiEfTr Si, Pz)0.^p^^\nf, Sf, - Qz)Ol^ , (15) 

where and 


Ot, — J dXF{ni, Si, X + Qt/ 2)'^^(}lF{nf, Sf, X — Qt/‘^) 

I KA [ ,1“ Sf \ 1 + ^2 

= 751 [A 1 

/ 1 lr>. 1 lf>_\l_V'n 

[loQo - i^Qz] 


, ,y I 1 + Si 

EM I Hi - -—,nf 

M {ni+ ^\ n/ - 


2)2 

1 + Sj \ 1 + 


... 1 + Sj 

+M ( Hi -- ,nf + 


2)2 

1 - Sf \ 1 - 


2 ' ^ 2 
In the above equation, M{nx,n‘ 2 ) is defined as 


Qy \ ( Qy 


l^Qy 


M{ni,n2) = J dxfn^ j fn^ “ 2 

= (2"i+"Vni!n2!)-'/2e-«?'/^ j dxe-^"Hr,, 

'Ql 


)± 

n2 



Qt 

‘71 

Qt 

71 


n2-n\ 2 

W'T 

e“ 4 r^ 2 -ni 

y 2 


ni-n2 -2 

e“ 4 irii-n2 

^n2 y 2 


Ql 


(ni < n2), 
(ni > n 2 ). 


(16) 


(17) 


where Hn{x) is the n-th Hermite polynomial, and L)^{x) is the associated Laguerre polyno¬ 
mial. 
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III. RESULTS 


In this section we show numerical results for pion emission from 1 — 2 GeV protons. In the 
numerical calculation we take the maximum Landau level to be ~ 3 x 10^, which somewhat 
limits our calculations. However, in order to test our model, we choose a proton energy of 1 
GeV and a strength of the magnetic held to be 5 x lO^^G. In this case the maximum Landau 
numbers are Urpm- = 45 for spin s = —1 and Umax = 50 for s = +1. Pion production has 


been shown 


22| to exceed photon emission in the region, y = = 0.01 ~ 1, in the 


semi-classical calculation. The value of y, 0.069, adopted in this calculation belongs in this 
region. 


In Fig. [T] we show the pionic decay 
widths of the proton at the Landau 
level rii and the spins state Si with 
a proton kinetic energy of 1 GeV. All 
results are summed over the hnal pro¬ 
ton spin and Landau levels. The solid 
and dot-dashed lines represent the de¬ 
cay widths of the proton with spin 
Sj = — 1 and Si = -|-1, respectively. 
For comparison we also plot results 
when the AMM is set to be zero, 
Kp = 0, with the dashed (sj = — 1) 
and dotted lines (sj = -|-1). 

When Kp = 0, the decay width of 
the proton with Sj = — 1 is a little 
larger than that with s* = +1, but 
these results are not signihcantly dif¬ 
ferent. In addition, the two states of 
s = ±1 are degenerate except for the 
lowest Landau state, and the spin- 
up and spin-down states cannot be 
uniquely determined when the AMM 


Decay Width of p to 



FIG. 1. (Color online) Decay width of a IGeV proton 
for the synchrotron emission of tt’s. The solid and dot- 
dashed lines represent the decay widths of protons with 
spin Sj = — 1 and Sj = 1, respectively. The dashed and 
dotted lines indicate results with s* = — 1 and s* = 1, 
respectively, for case that the anomalous magnetic mo¬ 
ment is omitted. 
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does not exist, Thus, this difference is 
not quantitatively meaningful. 

On the other hand, when the AMM is included, the width with s* = — 1 is enhanced up 
to about a factor of 100, while that with Sj = +1 is suppressed by a factor of about 1/100. 

In Fig|2l we compare our decay widths to those in the semi-classical approach. This figure 
shows the decay width averaged over the proton spin, T = \r{nmaxi +1) + ^{nmaxi ~l)]/2, 
when B = 5 X 10^®G as a function of the incident proton energy. 

When the proton energy is = 

1 GeV, the semi-classical approach gives 


Tsc ~ 180 MeV in Refs. 


MeV in 


22| and 210 


23|. These values are close to 


suits without the AIV 
expression of Ref. 


the maximum value of the pionic width of 
the proton with the AMM when Si = —1. 
However, the AMM is not taken into ac¬ 
count in the semi-classical calculations, 
and the results in the semi-classical ap¬ 
proach should be compared with our re- 



M. In addition, the 


23| is derived with the 


condition, y ^ 1, which is not consistent 
to the present condition, y 0.06 — 0.5. 
So, the decay width in the microscopic 
approach is much smaller than that in the 
semi-classical approach for the case with¬ 
out the AMM. Therefore, the AMM con¬ 
tribution, which is unique in the quantum 
mechanical approach, is vital for pion 
synchrotron radiation, particularly, in the 
limit of small curvature radii. 


Bp (GeV) 

FIG. 2. (Color online) Pionic decay width of protons 
versus the proton incident energy for B = 5x lO^^G. 
The decay width is averaged over the proton spin, 
and the Landau number is taken to be a maximum. 
The solid and dot-dashed lines represent the decay 
widths of the proton with and without the AMM, 
respectively. The dashed and dotted lines indicate 
the results in the semi-classical approaches of Refs. 


22l | and in 


23|, respectively. 


As the proton energy increases, the decay widths in the microscopic framework drastically 


increase, while those in the semi-classical approach increase more gently. In the semi-classical 
calculation one assumes that the Landau level number is very large. For example, in Ref. 23| , 













they take the limit, n* —?> oo. In Figj2l on the other hand, the maximum Landau is taken 
to be Umax ~ 20 — 700, which is too small for the semi-classical approximation. Thus, the 
semi-classical calculation is not justified based upon the present condition. However, the 
AMM contribution turns out to significantly increase the pion emission in the limit of a 
strong magnetic field. 



0.0 0.5 1.0 0.0 0.5 1.0 0.0 0.5 1.0 0.0 0.5 1.0 

COS0;r COS0;r COS0J COS0;r 

FIG. 3. (Color online) The differential proton pionic decay width versus the polar angle of pion 
emission. The widths are averaged over the initial Landau levels, 0 < Umax — ni < 4 (al,bl), 
5 < Umax — rii < 9 (a2, b2), 10 < Umax — ni < 19 (a3, b3) and 20 < Umax — Ui < 29 (a4, b4). The 
AMM is included in the upper panels (al—4), and not included in the lower panels (bl—4). The 
emitted pion energies are taken to be 100 MeV (dotted lines), 300 MeV (solid lines) and 500 MeV 
(dashed lines). 


In order to understand the angular distribution of emitted pions, in Fig. [3l we present the 
differential pionic decay widths of the proton, averaged over various initial Landau levels, 
0 ^ ‘^max ^2 ^ 4 (al,bl), 5 ^ '^max — 9 (^2, b2), 10 ^ ‘^max — 19 b3) 

and 20 < n^ax — Ui <29 (a4, b4). The AMM contribution is included in the upper panels 
(al—a4), but not included in the lower panels (bl—b4). The emitted pion energies are taken 
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to be 100 MeV (dotted lines), 300 MeV (solid lines), 500 MeV (dashed lines). Because the 
differential decay width is discrete for each angle, we average the results over the polar angle 
in every angular bin with A cos0^ = 0.1. 

We should note that the angu¬ 


lar distribution is symmetric between 
cosOt, < 0 and cos 6 '.,r > 0 , so that 
we do not plot results for all cos 6 *^. 
A pion produced from a proton with 
a higher Landau level Uj is emitted 
more transversely, and the AMM ef¬ 
fect shifts the pion emission to a more 
sideward direction when 0 < rimax — 
rii < 4. If we consider that the emit¬ 
ted pion decays into 27 s, it could af¬ 
fect the relativistic beaming of the 
energetic gamma rays from GRBs. 

Next we examine the proton-spin 
condition for the pionic decay width. 
In Fig. IHwe give the initial and hnal 
spin-dependence of the proton differ¬ 
ential pionic decay widths with (a) 
and without (b) the AMM. These 



- 1.0 


-0.5 


0.0 

COS0;, 


0.5 


1.0 


FIG. 4. (Color online) The differential pionic decay 
widths of protons with (a) and without (b) the AMM 
included. The widths are averaged over initial Landau 
numbers, 0 < Umax — n* < 9. The solid, dot-dashed, 
widths are averaged over the initial represents the results when 

Landau levels 0 < Umax — ni < 9. = ~^f = Si = Sf = 1 , and 

The solid lines represent the results ~ '®/ ~ respectively, where Si{f) indicates the 
when the initial spin s* = -1 and the initial (final) spin of the proton, 
hnal spin s/ = 1. The dot-dashed, 
dashed and dotted line indicate cases 

for which Sj = —Sf = 1, s* = s/ = 1, and Sj = sj = —1, respectively. 

When K,p = 0, the contributions from the spin-hip, s, = —Sf, are about 100 times larger 
than those of the spin non-hip, Si = Sf. Semi-classically the amplitude of the pion emission 
is proportional to < >• The pion emission along the 2 ;-direction is not allowed 
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by energy-momentum conservation, and < ippq(T±ipN > is the dominant contribution. Here, 
the spin-flip contributions become much larger than those from the spin non-flip reaction. 

When the AMM is included, only the contribution from Sj = —s/ = 1 is about 10,000 
times larger than those of the other channels, and the non spin-flip contributions are not 
much different between those with and without the AMM. 

As shown in Fig. [1] the AMM increases the decay width by 100 times for the case of 
Sj = 1 = —Sf = 1, and decreases it for the case of Si = —s/ = —1. When s* = —s/ = 1, the 
effects of the AMM and spin-flip are synchronized, and they enlarge the width by up to a 
factor of 10^. When Sj = —sy = 1, the two effects cancel and the width stays on the same 
order as that of the non spin-flip transition. 
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FIG. 5. (Color online) Pionic decay widths of protons (al, a2) and the pion transverse momentum 
(bl, b2) as functions of rij—nj—(sj—sy)/2. The AMM is included in the left panels (al, bl), and not 
included in the right panels (a2, b2). The initial Landau number is fixed to be nj -|- (1 — Si)/2 = 45. 


Next, we examine contributions from the final Landau level. In Fig. |5]we show the width 
(al and a2) and the pion transverse momentum (bl, b2) as functions of n, — ny — (sj — sj)/2, 
where the initial Landau number is fixed to be rij-l- (1 — Sj)/2 = 45. The results are calculated 
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with and without the AMM in the left (al, bl) and right panels (a2 and b2), respectively. 

Here, we note that the peak positions of the decay widths and the pion transverse mo¬ 
mentum are the same. When Kp = 0, the peak position is almost independent of the initial 
and final spin. When the AMM is included, however, the peak position of the non spin-flip 
transition is the same as that without the AMM, but it is shifted to smaller values when 
Si = —Sf = —1 and to larger values when Sj = —s/ = -|-1. 



FIG. 6. (Color online) The transition strength function M{ni,nf), Eq. ()17p . when n* = 45 -|- (1 -|- 
Si)/2. The solid circle and triangle show the peak positions in the cases Si = —Sf = —1 and 
Si = —Sf = 1, and the opened diamond indicates the case of the spin-flip Si = —Sf when Kp = 0. 

In order to study the AMM effect more clearly, we next examine the transition strength 
M{ni,nf) defined in Eq. f[T7)) . In Fig. [H] we show the n/-dependence of M{ni,nf) when 
TT-i = 45 -|- (1 -|- Si)/2. The solid circle and triangle in the inset represent the peak positions 
in the cases with s, = —Sf = —1 and s* = —Sf = 1 with the AMM included. The open 
diamond indicates the case of a spin-flip Si = —Sf without including AMM. 

M{ni,nf) shows an oscillating behavior, but only the strength after the last peak con¬ 
tributes to the results. In this region the strength rapidly decreases with increases oirii—nf. 
A pion cannot be produced in the free kinematics, B = 0, because of energy momentum con¬ 
servation. Under the influence of a magnetic held, momentum conservation is not satished, 
so that a pion can be produced in the kinematical condition far from the free kinematics. 
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where M{ni,n 2 ) rapidly decreases. 

The AMM gives a repulsive potential for Sj = — 1 and attractive for Sj = 1. The transition 
with Si = Sf = —1 introduces an additional energy to be consumed in the pion production. 
Thus, the difference between the initial and final Landau-levels is shifted toward a smaller 
number. In addition, the transition with Si = sf = 1 reduces the production energy. In this 
kinematical region a small difference between the initial and final Landau-levels signihcantly 
changes the transition strength. Therefore, the AMM plays the important roles of increasing 
greatly the pionic decay width when Si = Sf = —1 and to decrease it s* = s/ = 1. 



FIG. 7. (Color online) The dependence of the proton decay width on the pion energy for a magnetic 
field of 7? = 5 X 10^® G. The solid and dashed lines represent the results with and without the 
AMM. The proton incident energies are taken to be 1, 1.5 and 2 GeV, respectively. 

Next, we study the energy spectrum of the produced pions. We assume spherical symme¬ 
try in the momentum distribution of the initial protons and perform the angular integration 
of the differential decay width. In Fig. [7] we show the dependence of the proton decay width 
on the pion energy assuming a magnetic held B = 5 x lO^^G and the proton energies of 
cp = 1, 1.5, 2GeV, respectively. The solid and dashed lines represent the results with and 
without the AMM, respectively. 

As the proton energy increases, the decay width becomes larger, and the peak pion energy 
also increases. When the AMM is included, the peak height is about 59 times larger than 
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that when it is neglected for Cp = IGeV, and its ratios are 36 for Cp = l.SGeV and 30 for 
Cp = 2GeV; This difference between the peak height with the AMM and that without the 
AMM becomes smaller with increasing proton energy. 



0 200 400 600 

e.(MeV) 


FIG. 8. (Color online) The pion energy dependence of the proton decay width when the proton 
energy is Cp = IGeV. The solid and dashed lines represent results with and without the AMM. 
The magnetic fields are taken to be 1 x 10^®G, 2 x 10^®G, and 5 x 10^®G. Results increase with an 
increase of the magnetic field strengths. 

In Fig. [ 8 ] we show the dependence of the proton decay width on the pion energy for a 
proton energy Cp = IGeV with magnetic fields, B = 2 x 10^®G, 5 x 10^®G and 8 x 10^®G. 
The solid and dashed lines represent the results with and without the AMM, respectively. 

As the magnetic field increases, the decay width becomes larger, and the peak pion energy 
also increases. The peak height with the AMM included is about 93 times larger than that 
without the AMM for i? = 1 x 10^®G, and its ratios are 84 when B = 2 x 10^®G, 56 when 
B = 5 X 10^®G; the difference between the peak height with the AMM included and that 
without the AMM becomes larger as the magnetic field decreases. 

As noted above, a pion can be produced in conditions far from that of free kinematics. 
As the magnetic held decreases, the breaking of momentum conservation becomes larger, so 
that the effect of the AMM becomes more signihcant. 
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IV. SUMMARY 


In this work we have calculated the pion synchrotron radiation from high energy protons 
propagating in strong magnetic fields in a microscopic quantum field theoretical framework. 
We solved the Dirac equation in a strong magnetic field and obtained the proton propagator 
from its solution. Then, we derived the pionic decay width of the propagating proton in a 
fully relativistic and quantum mechanical way. 

Our results turn out to be compatible to those obtained by classical approaches. In 
particular, we hnd out that the anomalous magnetic moment has a very large effect which 
enlarges the emission rate by about 50 times, when the proton energy is IGeV and the 
magnetic field is 5 x 10^®G. In actual magnetars the surface magnetic field is known to be 
of order B ~ lO^^G . In the present method we did not perform a calculation for such 
a magnetic field strength because of the large number of Landau levels involved: a few 
thousand to a few million. 

As the magnetic field decreases, the AMM effect becomes larger. In a small magnetic 
field the decay width is very small, so that the proton energy must be larger to produce 
pions; if the semi-classical estimate is carried, the proton energy is expected to be at least 
about 100 GeV, and the maximum landau number should be a few hundred thousand. As 
the proton energy increases, on the other hand, the AMM effects diminish. Though it is not 
easy to estimate results for B ~ lO^^G, one expects the AMM effect to remain. 

As for future studies, since the emitted pion can decay into two gammas with some 
angular dependence with respect to the magnetic held, the secondary produced gamma ray 
may affect the photons from GRBs. More detailed studies are necessary for further discussion 
of that additional effect. Since the pions as well as the photons (including vector mesons) 
may be emitted in strong magnetic helds, the weak boson can also be produced from the 
propagating proton. This boson can decay into a pair of neutrinos. Detailed calculations 
regarding the neutrino pair production is also in progress. 


ACKNOWLEDGMENTS 

Work at NAOJ was supported in part by Grants-in-Aid for Scientihc Research of JSPS 
(26105517, 24340060). Work at the University of Notre Dame is supported by the U.S. 


15 


Department of Energy under Nuclear Theory Grant DE-FG02-95-ER40934. 


Appendix A: Proton Pionic Decay Width 


In this section we show the detailed expressions of the proton pionic decay width, Eq. m- 
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In the above equation Re is written as 
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Re can be written explicitly as 
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R(l, 1) = - {[AAf + s.Sfia, • a/)] [{Re • PfRiP^E + P]l) - P^P^l - {Mn + ^p/MR^Ql] 

+ SiSf{QE ■ Cli){QE ■ O/) [RjE — PfE + ‘d{MN + Kp/MpiY~\ 

+ {M^ + Kp/MRAfS, [{PR - PfE)iQoat - Qzk) - 2Ql{EiR - R,,a°)] 

+ (M^ + Kp/MRAsf [{PR - PkKQoa} - Q^aj) - 2Ql{Efa} - R^,aJ] } , 

fl(2,2) = 1 {[AtA, + s,Sf(a, ■ a/)| [(P,i • 
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+ SiSfiQi ■ ai){QL ■ ttf) [PIl - PfL + 2(MAr - Kp/M^f] 

-{Mr, - Kp/M^)AfSi [(P;i - P;j(Qoa^ - Q.a°) - 2gi(P,< - P,,a°)] 

-(M^ - Kp/M^)Asf [(P;i - PfL){Qoa} - Q^f) - 2Ql{Efa) - P/.aJ] } , 

n2 

P(3, 3) = { \PiL ■ PfL + — (kp/Mat)^] {AiAf + SiSfQi ■ aj) — SiSf{QL ■ cii){QL ■ o-f) 

-{kp/Mn) [Qo{AfSia- - AiSjUj) - Q^{AfSiOa° - AjS/aJ)] } . 
n2 

— ^Mi, {{Ei + Ef){AfSia^ — AiSfa^j) — {Piz + Pfz){AfSia^ — AjS/a^)} , 

P(4,4) = I [PiL ■ PfL + Mjf — (ftp/M^v)^] {AiAf + SiSftti ■ af) — SiSf{QL ■ ai){QL ■ «/) 

-{kp/Mn) [Qo{AfSiai - AjS/a}) - Qz{AfSia^i - As/a°)] } • 
n2 

+ —^M]\f {{Ei + Ef){AfSial — AiSfa^j) — {Piz + Pfz){AfSia^ — AiSfa^j){ , 

E{^, 2) = PixPfT {Ql + SiSfQi ■ Qf] — 2SiSf{QL ■ a,i){QL ■ 0-f)} , 

O'T 

E{^, 3) = {— [PiT(P/L ■ Ql) + PfT{PiL ■ Ql)] + SiSf{af ■ a,)] 

+ {PiT + PfT)SiSf{QL ■ Cif){QL ■ 0,i) 

+{PiT — PfT){Misf + Kp/Mr,) [AfSi{Qoai — Qza^) — AiSf{Qoa^f — Qza^f)] } , 

O'T 

PQ^ 4) = {— [PiT(P/L • Ql) + PfT{QL ■ Pil)] [A^/ + SiSfttf ■ Qi] 

+ {PiT + PfT)siSf{QL ■ af){QL ■ a-i) 

+{PiT — PfT){MN + Lip/M]\f) [AfSi{Qoai — QzO-^) — AiSf{Qoa^j — Qza^f)] } , 

P(2,3)=P(2,4) 

Ot 

= — {- [{PiL ■ QL)PfT + {Ql ■ P/I,)PiT] [AAf + SiSf{af ■ Qi)] 

+ {PiT + PfT)siSf{Q ■ af){Q ■ tti) 

+{PiT — PfT){—Mj, + Kp/Mr,)[^fSi{Qoai — g^a°) — AiSf [QoPf — Qza^f)]} , 
P(3, 4) = —QxPiTPfT {AAf + SiSftti ■ ttf) , 

where Pi(/)i, = (Pi(/), 0, 0, Pi(/>) and Ql = PiL - PfL = {Qo, 0, 0, Qz). 
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